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Metal whiskers often grow across leads of electric equipment causing short circuits and raising 
significant reliability issues. Their nature remains a mystery after several decades of research. It was 
observed that metal whiskers exhibit large amplitude movements under gentle air flow or, according 
to some testimonies, without obvious stimuli. Understanding the physics behind that movements 
would give additional insights into the nature of metal whiskers. Here, we quantitatively analyze 
possible mechanisms of the observed movements: (1) minute air currents, (2) Brownian motion due 
to random bombardments with the air molecules, (3) mechanically caused movements, such as (a) 
transmitted external vibrations, and (b) torque exerted due to material propagation along curved 
whiskers (the garden hose instability), (4) time dependent electric fields due to diffusion of ions, 
and (5) non-equilibrium electric fields making it possible for some whiskers to move. For all these 
mechanisms we provide numerical estimates. Our conclusion is that the observed movements are 
likely due to the air currents or electric recharging caused by external light or similar factors. 


I. INTRODUCTION 

Metal whiskers present a significant factor behind the 
reliability of many electronic devices. In particular, tin 
whiskers threaten many technologies using Sn based sol¬ 
der alloys. The mechanisms of metal whisker develop¬ 
ment and their properties remain insufficiently under¬ 
stood after almost 70 years of researchiU— Several pro¬ 
posed models of whisker nucleation and growth remain 
insufficiently quantitative and/or debatable. '~ 16 



FIG. 1: Snapshots of ~ 2.5 mm long tin whisker separated 
by time interval of ~ 0.05 s. The whisker is moved by a 
minute air current created by a gentle expiration towards the 
whisker. Courtesy of the NASA Electronic Parts and Pack¬ 
aging (NEPP) Program, 17 http://nepp.nasa.gov/whisker 

Here we discuss one of the long known yet poorly 
understood whisker properties: their ability to exer¬ 
cise movements of significant amplitudes comparable to 
whisker lengths. These movements have been observed 
multiple times in response to rather gentle (say, expira¬ 
tory) air flows^ 7 ’ 1 ® as illustrated in Fig. [L] 

More intriguing is that many researchers have re¬ 


ported spontaneous (unintended) whisker movements 
taking place without any obvious stimuli. The author too 
has observed unintended movements of some Zn whiskers 
using a magnifying glass and a LED flash light (whisker 
infested Zn samples, courtesy of J. Brusse, the NASA 
NEPP Program); however the question of reliability of 
those observations or possible uncontrolled air flow ef¬ 
fects remains open. 

The testimonies of ‘spontaneous’ metal whisker move¬ 
ments have been a topic for extensive discussions in the 
metal whisker community, including through the ‘tin 
whisker group teleconference’^ which resulted in vari¬ 
ous hypotheses, such as, e. g. ‘Brownian movements’, 
‘mechanical vibrations’, or ‘garden hose instability’, and 
some others. These hypotheses about metal whisker 
movements, either air current stimulated or spontaneous, 
have not been quantitatively tested. 

This paper provides a so far missing quantitative anal¬ 
ysis of whisker movement mechanisms. Its results be¬ 
low are accurate to the order of magnitude. This ap¬ 
proximate approach appears adequate in view of various 
uncertainties related to whisker parameters, as well as 
the state of the art approximate nature of theoretical 
descriptions of stochastic and hydrodynamic phenomena 
involved. The latter issues make exact model solutions 
excessive and beyond any reasonable accuracy. 

Assuming a diverse readership including many non¬ 
physicists, a comment may be in order regarding the na¬ 
ture of the order-of-magnitude approximate analysis. It 
is well known with many subjects in the hydrodynamic 
theory, statistical physics, the theory of disordered sys¬ 
tems, and many other subjects where exact solutions 
are impossible or useless. An important lesson learned 
from that multitude of applications, is that the order-of- 
magnitude approach provided results that are generally 
good to the accuracy of insignificant numerical factors, 
such as, say, 0.7, or 2.1, etc., as compared to the exact 
results where available. In the meantime, it is more in¬ 
tuitive, fast, and economical avoiding unnecessary math¬ 
ematical complications. 
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Following the standard notations, the equality sign = 
will be everywhere replaced with the ‘equal in-the-order- 
of magnitude’ sign ~, meaning ‘to the accuracy of a nu¬ 
merical factor’. For example, the area of a circle of di¬ 
ameter d is estimated as A ~ d 2 , instead of the exact 
A = nd 2 /4, different by a numerical factor of 7r/4 ss 0.75. 
Below, we will compare, at some points, the ‘order-of- 
magnitude’ and exact results. 

In addition, the ‘order-of-magnitude’ analyzes often 
use intuitive estimates for derivatives and integrals, such 
as dy/dx ~ Ay/Ax and J q “ +A ' Y y{x)dx ~ y(a) Ax where 
Ay and Ax are the intervals for variables x and y. These 
approximations generate additional numerical uncertain¬ 
ties. For example, assuming y = x 2 , the latter estimates 
yield, dy/dx ~ x 2 /x = x and ydx ~ yx = a; 3 , missing 
respectively the numerical factors of 2 and 1/3, insignifi¬ 
cant as not changing the order of magnitude of the esti¬ 
mated quantities. Many useful examples of such approx¬ 
imations are described in the classical text by Migdah— 

The observed insignificance of numerical factors, typ¬ 
ically not too different from unity, has found a com¬ 
monly shared conceptual explanation (folklore attributed 
to Einstein) that, because numerical uncertainties at dif¬ 
ferent steps in complex processes are mutually indepen¬ 
dent, they tend to cancel each other similar to random 
displacements in diffusion processes (i. e. ‘errors prop¬ 
agate by diffusion’). Say, the coefficients of 2 and 1/3 
in the above examples, can combine into 2/3, which is 
even closer to unity. There are very few known excep¬ 
tions where such cancelation does not take place. This 
amazing observation is fully consistent with the fact that 
the laws of physics typically do not include any signifi¬ 
cant numerical multipliers when written a ‘natural’, say, 
Gaussian units (as opposed to SI where such multipliers 
are enforced). 

Using the ‘order of magnitude’analysis throughout 
multiple subjects in this paper will offer a practicum that 
can be helpful to graduate students and other audience 
interested in practical approximate techniques of theo¬ 
retical physics. In addition, this paper provides a brief 
tutorial on general bending dynamic properties of elas¬ 
tic rods (modeling metal whiskers); they are not often 
taught in the standard physics curriculum and can be 
useful for researchers interested in metal whisker effects. 

The paper deals with multiple independent subjects; 
hence, quite a few quantities denoted with various letters. 
The intent was to keep the notations close to standard, 
such as, e. g. A for area, D for diffusion coefficient, etc. 
One obviously inconvenient case is d for the diameter, 
which meaning could be mixed with that of differential. 
To avoid any confusion, d is never used as a differential in 
what follows; in several cases where the notion of differ¬ 
ential becomes unavoidable, we use the partial derivative 
denoted with d. 

The paper is organized as follows. Section HD describes 
the elastic beam model of a metal whisker, and espe¬ 
cially its bending properties for uniform and nonuniform 
whiskers. Sec. Mil provides estimates of possible whisker 


movements due to viscous air flow. In Sec. hyi the anal¬ 
ysis of Brownian movements is done for metal whiskers 
in comparison with the standard Brownian movement of 
small particles, as well as the viscous drag effects. Pos¬ 
sible mechanical causes of whisker movements are ana¬ 
lyzed in Sec. |V] Finally, Sec. |Vj discusses the possibility 
of whisker movements due to electrostatic effects. The 
numerical estimates corresponding to all the analyzed 
mechanisms are summarized in Sec. IVIII The paper con¬ 
clusions are presented in Sec. IVIII I 


II. ELASTIC BEAM WHISKER MODEL 

Because the observed whisker movements are re¬ 
versible, the elastic approximation appears to be ade¬ 
quate. We describe each whisker as a cylinder shaped 
elastic beam of a relatively small diameter, d <C l, where 
l is the whisker length, parallel to the x-axis of a Carte¬ 
sian coordinate system (Fig. [2]) (a)), and limit ourselves 
to bending deformations. 


A. Euler-Bernoulli equation 


The elastic beam deformations are described in terms 
of deflection u(x,t) by the EulerBernoulli equation (see 
e. g. Refs. I21II22I or dozens of other excellent sources) 


^ rT d A u d 2 u du 

YI w +r ’ A W + B Tt =,{x ’ t) 


(i) 


where Y is the Young’s modulus, I is the polar area mo¬ 
ment of inertia of the cross-section (defined as the inte¬ 
gral of the square of a radius vector of a typical point in 
the cross-section over the cross-section area), p is mass 
per unit volume constant throughout the cross-section, 
A is the cross-sectional area, and B is the friction force 
coefficient, and f{x,t) is the distributed force per length 
that depends on the coordinate and time. 


B. Qualitative description 

While Eq. (JT]) should describe the entire beam dynam¬ 
ics, in what follows, we rely on an approximate ‘order-of- 
magnitude’ approach facilitating its understanding, and 
neglecting possible quantitative discrepancies in numeri¬ 
cal coefficients. 


1. Uniform whiskers 

Here, we consider whiskers with uniform cross-section 
along the longitudinal coordinate. We start with estimat¬ 
ing the beam deflection A in response to force F, either 
distributed or exerted close to the beam edge. Consider 
an ark of length l illustrated in Fig. [2] (b). Its upper 
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FIG. 2: (a) Elastic beam of length l and diameter d char¬ 
acterized by the deflection u(x, t) from axis x. (b) A beam 
bent with the curvature radius R and maximum deflection 
f)]max = A (not to scale d -C l). Fat tangential arrows 
show stresses that have opposite directions above and below 
the neutral line presented with a dash dot line. 


and lower parts are respectively stretched and squeezed 
due to bending. A line between them, called the neu¬ 
tral axis, is neither stretched no squeezed maintaining 
the same length as before deformation. Denoting the arc 
curvature radius by R and using the geometry in Fig. [2j 
it is straightforward to estimate the maximum deflection 
as 


A - l 2 /R. (2) 

The length of arc at distance d above (below) the neu¬ 
tral line is (R ± d)a where a is the arc angle. There¬ 
fore, the relative deformations of the upper and lower 
parts are estimated as ±A l/l ~ ±d/R. The correspond¬ 
ing stresses are s = YAl/l ~ ±Yd/R where Y is the 
Young’s modulus. They result in the torque density ±sd 
(i. e. force density s times lever d), and the integral 
torque M ~ Asd ~ Yd 4 /R, where A ~ d 2 is the cross 
sectional area. Expressing R from Eq. ©, one gets the 
torque related to the bending deformation, 

M~Y^-. (3) 

On the other hand, the torque can be expressed 
through a more or less perpendicular exerted force F as 
M ~ FI. Comparing the latter to Eq. © yields, 


As an illustration of the ‘order-of-magnitude’ ap¬ 
proach, the exact result^ is that force F acting per¬ 
pendicularly on the end of a cantilever beam creates a 
displacement given by Eq. © with a numerical coeffi¬ 
cient of 1/3, while the case of uniformly distributed force 
yields the coefficient 1/6. 

Eq. © shows that bending deformation A is by the 
factor ( l/d ) 2 1 greater than the compressive deforma¬ 
tion F/K. In other words, the bending elastic spring Kb 
is by the factor ( d/l ) 2 -C 1 smaller than the compressive 
spring constant K. 

The above simple estimates lead to the characteristic 
frequency of transverse vibrations, 



where m ~ pld 2 is the beam mass. We observe that 
bending vibrations have frequencies by the aspect ratio 
d/l <€. 1 lower than the standard elastic vibration fre¬ 
quencies LO. 

Estimates in Eq. © can be verified in several ways. 
For example, by setting equal the terms on the left-hand- 
side of Eq. © and estimating d 4 u/dx 4 ~ u/x 4 and 
d 2 u/dt 2 ~ uu) 2 . Or by considering the standard angular 
dynamics via Iod 2 (j)/dt 2 = M where 

Io ~ ml 2 ~ pl 3 d 2 (6) 

is the moment of inertia, 4> is the deflection angle, and 
M is given by Eq. © with A ~ l<j>. 

To further illustrate the accuracy of the ‘order-of- 
magnitude’ treatment, note that the vibrational spectra 
of cantilever beams are described as^i 


u> b 




( 7 ) 


where ft 2 is a numerical factor labeling different vibra¬ 
tional modes (in terms of the number of nodes). Taking 
into account that / ~ A 2 , Eq. © coincides with Eq. © 
to a numerical factor /3 2 , which for the first mode is close 
to one. The coefficient /3 increases with the mode number 
not accounted for by the approximate Eq. ©. However, 
the existence of higher excited vibrational modes is irrel¬ 
evant to this paper subject not dealing with resonance 
phenomena. 

Since 


FI 3 F / l\ 2 _ F 
Yd 4 ~K\d) ^~K h 


( 4 ) 



( 8 ) 


Here, we have introduced the compressive spring con¬ 
stant K = Yd, 2 /l as the proportionality coefficient be¬ 
tween the force sA and compressive elastic deformation 
A l, and the effective bending spring constant Kb = 
K(d/l) 2 . 


represents the sound velocity v s , one can conclude that to 
is the characteristic reciprocal time, over which the sys¬ 
tem retains a coherent translational dynamics (all parts 
of the body move coherently). u>b represents a similar 
value for the bending dynamics. 
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The latter can be shown more explicitly by approxi¬ 
mating F in Eq. © through the Newton’s second law 
as 


F ~ roA/r t 2 with m ~ pcPl (9) 

where tj, is the deflection time and the acceleration is 
estimated as A/r 2 . Substituting Eq. © into Eq. © 
yields, 


n ~ 


d_ 

v s 



1 

U>b ' 


( 10 ) 



The estimates in Eqs. © and © are limited to dy¬ 
namics governed by the whisker bare mass m. They do 
not apply to the case of viscous drag driven whiskers, par¬ 
ticularly in Brownian motion, where the dynamic prop¬ 
erties are determined by a composite system of a whisker 
and a fluid (see Sec. hyi) 

Consider next the viscous drag on a cylinder whisker 
with high aspect ratio l/d 1 that corresponds to the 
third term on the left hand side in Eq. ©■ Because d 
is the only characteristic length at hands, one concludes 
that the air velocity field v a (r) is significantly disturbed 
at distances r < d from the whisker; hence, |Vu a | ~ 
v/d where v is the whisker velocity. As a result, to the 
logarithmic accuracy*^ 24 the viscous drag is estimated 
as 


F ~ r]A s \S7va\ ~ rjvl (11) 

where A s ~ Id is the whisker side area. 

For completeness, we briefly describe also the condi¬ 
tion of buckling instability. It represents switching from 



FIG. 3: Side and top views of an airflow around a long cylinder 
shaped whisker with arrows showing the current lines. F is 
the integral force on the cylinder. 


FIG. 4: Buckling instability results in switching from verti¬ 
cal beam purely compressed under force Fo (left) to a bent 
configuration characterized by a finite deflection A (right). 


a purely compressed (no beam deflection) to a laterally- 
deformed state characterized by a finite deflection A il¬ 
lustrated in Fig. U The switching between these states 
takes place when the compressive force Fo exceeds a cer¬ 
tain critical value Fo c that depends on the beam length 
l, or, in other words, when, the beam length l exceeds a 
certain critical values l c that depends on Fo. 

To derive the switching condition, we note that the 
bending force F that appears in Eq. © can be expressed 
through the originally applied compressive force Fo as 
F = Fo sin cf> where <f> is the bending angle shown in Fig. 
21 Estimating (j) ~ X/l and substituting F = F 0 \/l in 
Eq. © yields, 

F 0 = F 0c ~ Yd 4 /l 2 . (12) 

When F 0 < F 0c , Eq. © cannot be satisfied. On the 
other hand, when Fo > Fo c , the above estimate for the 
angle becomes inaccurate and we are left with Eq. ©. 

Expressing the beam length through Eq. m , gives 
the critical beam length, 

lc = s/Yd 4 /F 0 . (13) 

The above described buckling yields a significant deflec¬ 
tion A < (. The results in Eqs. m and m coincide 
with that of exact analysi s 21 ' 22 to the accuracy of numer¬ 
ical coefficients of the order of unity. 


2. Nonuniform whiskers 

Consider next whiskers with significantly nonuniform 
cross-sections i. e. with diameters d(x) significantly vary¬ 
ing along the whisker longitudinal coordinate x, as illus¬ 
trated in Fig. [5] Such a whisker can be thought of as 
a set of small domains, such that d(x) in each of them 
is more or less uniform. In static equilibrium, the forces 
on subsequent domains must be equal. Therefore, Eq. 
© predicts then that the bending angles in individual 
domains are inversely proportional to d 4 , i. e. a whisker 
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integral bending is determined by its constriction region 
(if any) where the diameter d c -C d. 

More specifically, taking the individual domain lengths 
Sli = di = d(xi ), setting subsequent domain forces equal 
and using Eq. © yields 


59 , 


K_ J_ 

5k K dl 


(14) 


Eq. (Tl4l) can be further generalized to give the integral 
bending angle, 


(15) 

i 1 

where p(di) is the probability to find a domain of diam¬ 
eter di, so that YliP{di) = 1- 

As an example, consider a whisker of length l = 1 
mm and average diameter d av = 10 //in having a con¬ 
striction with diameter d c = 1 //m extending over dis¬ 
tance 51 ~ d av corresponding to the sketch in Fig. [5] 
In the constriction region, we set 5k = 1 /tm which 
makes n c ~ 10 domains forming the constriction and 
yielding p(d c ) ~ d av /l ~ 0.01. On the other hand, the 
average diameter region will contain n av ~ l/d av = 10 2 
domains of length 5k = d = 10 //in each. According 
to Eq. (lisp , each domain in the constricted region will 
contribute the deflection angle 0i that is by the factor 
(d av /d c ) 2 = 100 greater than that by the average diam¬ 
eter domain. Therefore, the relatively short constriction 
region (of ~ 1 % of the total length) will contribute to 
the integral bending, approximately the same angle as 
the rest 99 % of the whisker, as dictated by the relation 
n c /d 2 c « n av /d av . 

It is obvious that constrictions will affect whisker vi¬ 
brations. This effect can be readily described for the 
case of a distinct local constriction of the type presented 
in Fig. [5] The constriction related bending elasticity, 
Kb c is softer than that of the average one, 

l d 2 Id 2 

Kbc ~ y%t K b ~ ytK h « K b , (16) 

'•c w av u av 


* 



FIG. 5: Sketch of a whisker with constriction where the di¬ 
ameter d c is much lower than its average diameter d. 


while the vibrating mass m is more or less the same. 
As a result, the characteristic vibrational frequency cor¬ 
responding to a vibrating whisker with constriction is 
estimated as 



(17) 


and the bending time r bc is by the same factor longer 
than that (17,) without the constriction. 


III. AIR FLOW EFFECTS 


Here we test the hypothesis of minute current flows 
causing the observed whisker movements. The fact that 
whisker deflections A can be significant, is readily quan¬ 
tified by using Eqs. © and m- Combining them and 
setting A ~ l results in the following estimate for the lo¬ 
cal air velocity that produces significant (comparable to 
the whisker length) deflections, 


v 


Yd 4 

rjl 3 


(18) 


The numerical estimates in Sec. IVIII show that thin 
(micron in diameter, or constricted to such neck diame¬ 
ter) and simultaneously tall (l 1cm) whiskers can be 
significantly moved by rather modest air currents of the 
order of 1 mm/s. On the other hand, thicker (d ~ 10 //m) 
whiskers require unrealistically strong winds v a > 10 4 
cm/s to be moved. 

Table Q] below shows for comparison some realistic fig¬ 
ures for the air flow velocities in the typical lab settings. 
They all are sufficient to create movements of thin (or 
constricted) metal whiskers. This is further illustrated 
in Table Hill below. 


TABLE I: The characteristic air flow velocities due to various 
sources. 


Air flow cause 

Air flow velocity (cm/s)“ 

Room A/C 

400 

HVAC Vent 

220 

Person walking 

180 

Door opening 

120 

Expiratory air flow 6 

100 

Diffuser vent 

25 


“Data from Ref. [25I . except the expiratory mode, for which 
the velocity was estimated assuming the average lungs expiratory 
volume^ of ~ 1 L and the lips opening area ~ 1 cm 2 held for 10 s. 

b Also mentioned in the NASA whisker databasa2i as the “gentle 
air flow” produced by pursing one’s lips and gently blowing in the 
direction of the whisker 
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IV. BROWNIAN MOVEMENTS 

Here, we consider the Brownian movements of whiskers 
starting with the classical analysis for a suspended par¬ 
ticle in a fluid or gas of viscosity p. The rationale for 
this consideration is that the reported chaotic whisker 
movements seem visually similar to the Brownian pat¬ 
terns, and that whisker diameters are in the range of the 
typical Brownian particle dimensions. 


A. Classical Brownian movements 


The analysis is based on the characteristic relaxation 
time Tjj of a particle of diameter d and material density 
p ~ m/d 3 . The viscous drag force F ~ pvd alters the 
particle momentum by mv over time t v . From that, one 
estimates, 


m 




(19) 


T rj escribes the time interval during which a particle main¬ 
tains a coherent translational motion. 

A particle trajectory is then considered as consisting of 
mutually independent steps of duration t v each. A one 
step displacement, 


a ~ vtt, 7 



y/k B Tpd 

V 


( 20 ) 


constraints In particular, for the case of a harmonic 
oscillator, the classical result L oc y/t holds for relatively 
short displacements L < Lt ~ yJk B T /K , where Lt is 
the time independent equipartition value limited by the 
restoring force F = —KL. 


B. Brownian movements of metal whiskers 

Along the same lines, one can consider random whisker 
bending over angle <f> during time t as a succession of 
t/T v 1 events, each producing the characteristic bend¬ 
ing angle A <f> ~ lott v where lot is the thermal angular 
velocity and r r/ is the coherency time. 

Unlike the classical Brownian motion, the stochastic 
displacements here are limited by the elastic bond setting 
the maximum ‘equipartition’ value of whisker deflection 
A t- This limitation is similar to that of Brownian move¬ 
ments of a harmonic oscillator mentioned at the end of 
Sec. II V Al Using the equipartition theorem in the form 
of KbXj. ~ fcsT, yields the deflection At and its corre¬ 
sponding angle A(f> = At/ l, 


A cj) ~ 



( 22 ) 


The characteristic thermal angular velocity correspond¬ 
ing to the equipartition theorem lot ~ y/k B T/I q is self- 
consistently related to A</>, 


where we have used the characteristic thermal velocity 
vt ~ \Jk B T/m, k B is the Boltzmann’s constant, and T 
is the temperature. During time t r^, the particle 
makes t/T v 1 random steps resulting in the average 
displacement 





lk B T 

pd 


t = yj D d t with D d 


k B T 

pd 

( 21 ) 


Here D d is the particle diffusion coefficient. 

To the accuracy of numerical multipliers, Eq. 
m coincides with the known results for Brownian 
movements ,S 9 ~ 3 2 As yet another illustration of this pa¬ 
per ‘order-of-magnitude’ approach, note that the exact 
result for the Brownian movement amplitude is different 
by the numerical coefficient ~ 0.5 from that in Eq. m- 
The characteristic amplitudes of Brownian motion are 
rather small, below the naked eye range. Indeed, us¬ 
ing the viscosities of air and water to be respectively 
p ~ 10 -4 g/cm-s and p ~ 10 -2 g/cm-s, yields respec¬ 
tively t v ~ 10~ 3 s and t v ~ 10 -5 s, and the diffusion 
coefficients D d ~ 10~ 6 cm 2 /s and D d ~ 10 -4 cm 2 /s. On 
the typical observation time scale f ~ Is, for a micron 
sized particles, d ~ 1 /rm, the characteristic displace¬ 
ments will be respectively L ~ 1 /iin and L ~ 10 /im, in 
agreement with the results of optical microscopy^ 

The concept of classical Brownian movements is sig¬ 
nificantly modified by the presence of external fields or 


A0~wr/wfc. (23) 

Note that Eq. (1221) does not imply that the fluctu¬ 
ation A <f> is produced by a ‘rigid body’ movement of a 
whisker as a whole. The same result can be obtained 
by considering a whisker made of mutually independent 
small domains of length 51 <C l each. A single do¬ 
main bending dispersion is similar to that in Eq. (1221) . 
(< 5(j)) 2 = ( kT/Yd 3 )6l. Multiplying ( 5q i) 2 by the number 
of such domains l /51 gives the dispersion A</> that repro¬ 
duces Eq. (l22l) . 

In connection with the latter paragraph, note that the 
characteristic angular velocities of individual domains are 
by the factor ( 1/51 ) 3//2 1 greater than the integral an¬ 
gular velocity lot- Hence, the latter derivation implies 
high frequency chaotic dynamics with rapidly varying 
local configurations limited to the integral angle of Eq. 
( 1221 ) . In fact, Eq. ( 1221 ) reproduces the classical result for 
fluctuations in curvature of long molecules^ 

We now estimate the characteristic time of establishing 
the above maximum deflection. Similar to the treatment 
of classical Brownian movements in Sec. IIV Al the angle 
of deflection over time t is accumulated via t/r^ mutually 
independent characteristic partial deflections 5<j) ~ lott v , 
i. e. 


A 



(24) 
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To evaluate r r/ , we use the Newton’s second law pre¬ 
senting the viscous drag of Eq. m in the form F ~ 
mv/T v . This yields t v ~ pd 2 /rj, which coincides with the 
estimate for a sphere of diameter d in Eq. (flUl) (again, 
to numerical factors). 

Combining Eqs. E21) and flU and substituting t v 
yields the characteristic time of reaching the maximum 
bending Aej) [of Eq. 

< 25 > 

Numerical estimates in Table mu below show that, for 
realistic parameters, the predicted Brownian bending an¬ 
gle Aej) is too small to explain the observed whisker move¬ 
ments, although possible bending times of the or¬ 
der of seconds are compatible with the data. Even for 
severely constricted whiskers with local diameter d ~ 0.1 
/iin, the predicted bending angle Aej> oc d~ 2 remains small 
compared to the observations. 


C. Brownian movements vs. viscous drag 


The two mechanisms are additionally illustrated in Fig. 
[6] where a particle is bombarded with molecules within a 
layer of thickness Iq . We consider the particle movements 
in the vertical direction, parallel to which, there is a flow 
of air with velocity v a at distance ~ d in horizonal direc¬ 
tion (d is the characteristic distance over which the air 
flow is disturbed by the particle). The airflow velocity 
at distance Iq from the particle is estimated as ~ v a lo/d. 
Each molecule of mass mo hitting the particle in hori¬ 
zontal direction, transfers the momentum moV a lo/d, and 
the momenta transferred from the left and right do not 
cancel each other. The momentum transferred over time 
t is estimated as 

Ap visc ~ nv T Atm 0 v a l 0 /d 

where A is the particle area. Dividing the latter by time 
t results in the viscous drag force rjv a d with the stan¬ 
dard viscosity r] = nloksT/vr • To the contrary, for the 
Brownian movements, the momenta transferred from the 
top and bottom collisions mostly cancel each other and 
only the 1/y/N fraction of them remains unbalanced as 
a fluctuation (we skip the details of derivation). 


A comment is in order regarding one intuitively jus¬ 
tified hypothesis, that spontaneous whisker movements 
could be of the same nature as that of the often observed 
stumbling of small dust/smoke particles in the air^ both 
being Brownian. Here we briefly discuss how Brownian 
movements are typically much weaker than that due to 
viscous drag, and are hardly relevant for metal whiskers. 

The ratio of the characteristic displacements due to 
viscous drag, v a t over that due to Brownian diffusion, 
V Ddt can be expressed in terms of microscopic parame¬ 
ters if we take into account the definition for Dj in Eq. 
m along with the standard representation for viscos¬ 
ity, p ~ hUbTIq/vt where vt = y/ksT/mo is the ther¬ 
mal velocity, mo is the molecular mass, n is the molecule 
concentration, and Iq = 1 /{no) is the molecule mean free 
path with a being the molecular scattering cross section. 
This yields, 


Vat ^ Va_ VTtd 

y/D d t vt V cr 


(26) 



FIG. 6: A discrete (4-side) sketch of a particle bombarded 
by molecules having on average the thermal velocity vt and 
capable of colliding with the particle when they are within 
the range of the mean free path lo from it. Shown in short 
arrows are vertical (in the plane of figure) air currents of low 
velocity v a that is reduced by the factor lo/d in the domain 
of molecules capable of colliding with the particle. 


The right-hand-side in Eq. (E6l) is readily estimated. 
We take a modest v a ~ 1 mrn/s, the typical vt ~ 10 5 
cm/s and cr ~ 10~ 16 cm 2 , also setting d ~ 1 pm, and 
i ~ Is corresponding to the naked eye observation times. 
As a result, one gets, v a t/y/Ddt ~ 300. The viscous 
drag will be yet more significant for bigger particles and 
stronger air currents. 

It may be worth pointing at a more intuitive interpre¬ 
tation of the above. Presenting 1/cr = Ion translates the 
square root in Eq. (l26l) into the form y/~N = y/nlovtd 
with N 1. Therefore, this is a very large number 
N 1 of molecules involved that makes the viscous drag 
more efficient by the factor proportional to y/~N. 


Our conclusion here is that the Brownian motion is 
hardly observable with metal whiskers. Similarly, the 
most visible movements of dust particles in the air, are 
often explained by tiny air currents due to thermal con¬ 
vection or other ill controlled factors. 


V. MECHANICAL MODELS 
A. Sample vibrations 

Here we test the hypothesis that metal whisker move¬ 
ments are caused by uncontrolled vibrations of the film 

































transmitted from remote sources, such as some power 
equipment, transportation, etc. 

The characteristic deflection caused by such vibrations 
can be estimated by using the non-inertial frame of ref¬ 
erence pinned to the sample, as illustrated in Fig. E 
In that frame, there will be the inertial (fictitious) force 
F t = —mw exerted on the whisker, where w is the sample 
acceleration relative to the laboratory system. Substitut¬ 
ing that force in Eq. a and using m ~ pcPl, yields 


pd 2 w (l \ 

w 


(27) 





FIG. 7: Left: whisker movements due to the inertial (ficti¬ 
tious) forces caused by the sample mechanical vibrations with 
acceleration w. Dash-dot arrow illustrated the sample accel¬ 
eration. Right: whisker deflection caused by the momentum 
change due to its growth (garden hose model). The dashed 
arrows represent the material movement to the base of the 
whisker and its caused movement (growth) of a whisker as a 
whole. 

For numerical estimates in Table cm we used the 
extremum value related to acceleration due to gravity, 
w = g = 9.8 m/s 2 (that can be felt in a roughly mov¬ 
ing truck—). Both that and a ‘slightly rough— accel¬ 
eration w ~ 0.1 < 7 , are well above the standard lab ex¬ 
posed ‘smooth’ accelerations, which fall in the ranged 
of w < 0 . 01 < 7 . For the latter, the whisker deflection are 
at least by the factor of 0.01 smaller than shown in Ta¬ 
ble IIIII (for w = g), and they are much smaller than 
observed. However, severely constricted whiskers with 
d ~ 0.1 pm can respond by deflections that are by two 
orders of magnitude higher than given in Table ITOl They 
would be consistent with the observed movements. 

A word of caution is in order regarding possible reso¬ 
nance phenomena amplifying the vibrational amplitudes 
for frequencies in the observed whisker motion range, say 
cob ~ 30 in Table IIIII While the possibility of such ampli¬ 
fication depends on the spectral composition of a noise, 
the resonance amplification factor for whiskers can be sig¬ 
nificant, uiT r 1 where the velocity relaxation time due 


FIG. 8 : Charge patch model. Left: positive and negative 
charge patches of characteristic linear dimension do generat¬ 
ing random electric field vectors at distances l 3> do- Right: 
electric field fluctuations with absolute value decaying as 
E oc l/l with distance l do from the surface. 


to viscous friction is estimated as r r ~ m/(j]d). Using the 
numerical parameters from TablclIIlvields cot,, ~ 10—100. 


B. Garden hose instability 

The fluid-conveying pipe, particularly a garden hose 
can exhibit oscillating behavior due to interactions be¬ 
tween the moving fluid and a hose. Various mechanisms 
of such interactions have been considered (see Refs. I3l39l 
and references therein). Here we concentrate on one of 
them applicable to the very slow material flow relevant 
for the case of growing whiskers. 

Denoting v w = dl/dt the time independent characteris¬ 
tic whisker growth, the corresponding momentum trans¬ 
fer is due to the change in whisker mass, 

dp dm l9 dl 9 ,9 . . 

M =v "-M~ v - pd ai = l ’ v '- d - (28) 

Assuming the whisker bending angle (~ 90°) such that 
the entire force F = dp/dt is directed more or less per¬ 
pendicular to the whisker, and substituting that force in 
Eq. gD yields the relative whisker deflection, 



The latter ratio is extremely small because of the rela¬ 
tively low whisker growth velocity ~ 1 A/s compared 
to the sound velocity v s ~ 10 5 cm/s. As seen from Table 
mu this mechanism should be totally rejected even for 
severely constricted whiskers. 
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VI. ELECTRIC FORCES 

Here we consider a hypothesis that whisker move¬ 
ments are caused by interactions of whiskers with charged 
patches on the underlying surface. 

A. Charge patch model and static deflection 


needed. Here we follow the concept of electric field fluc¬ 
tuations due to the time dependent electric dipoles on 
metal surfaces (or in its underlying layers) . 46 ' 47 

The dipoles are caused by ions diffusing across the sur¬ 
face. Each diffusing ion of charge e creates a mean square 
dipole fluctuation ( Sp) 2 ~ e 2 Dit where Di is the ion diffu¬ 
sion coefficient. The corresponding dispersion in electric 
fields at distance l from the surface becomes 


The charge patch model^ takes into account local 
electric charges in the metal surface originating from var¬ 
ious imperfections: grain boundaries, local deformations, 
variations in chemical composition, contaminations, dis¬ 
locations, etc. It was able to describe many properties 
of whiskers, such as the existence of dormant and lin¬ 
ear growth stages, their high aspect ratios, growth ve¬ 
locities, versatility of factors affecting whisker appear¬ 
ances and kinetics, and broad statistical distributions of 
whisker parameters. That model however calls upon fur¬ 
ther verifications. Earlier evidence of electric bias ef¬ 
fects on whiskers did not discriminate between the field 
and current effects.—^ However, movements of metal 
whiskers caused by the electric field was experimentally 
demonstrated^, and there is a preliminary evidence that 
the electric field strongly accelerates whisker growth. 43,44 

The charge patch model assumes relatively small sur¬ 
face regions illustrated in Fig. [5]with characteristic linear 
dimension do 10 A™ where individual patch charges are 
mutually uncorrelated. In the near surface region l < do, 
the surface field is more or less uniform and perpendic¬ 
ular to the surface; its absolute value Eq is determined 
by the local charge density. Farther from the surface, at 
l ^ do, the oppositely charged random patches mostly 
cancel each other contributions making the field decay as 
E ~ Eodo/l on average. In that latter region, the field 
vectors are oriented almost randomly. The characteristic 
maximum surface field at l < do was estimated to fall in 
the range of Eq ~ 10 4 — 10® V/cmji^ 

For a whisker of length l at angle 6 with the normal, 
the characteristic force parallel to the surface can be esti¬ 
mated as F ~ (pE/l) cosO ~ (pE 0 do/l 2 ) cosd where the 
whisker dipole moment is p = (3E and its polarizability 
/3 is estimated as that of a long thin metal rod^ i. e. 
/3 ~ l 3 . Substituting that force in Eq. 0 yields 



Numerically, the latter ratio can be comparable to (or 
even formally exceed) unity for reasonable parameter val¬ 
ues as illustrated in Sec. IVIII Therefore, the electrostatic 
forces due to charge patches are strong enough to cause 
a significant static whisker bending with A /l < 1. 


B. Ion dipoles 

To estimate the possibility of whisker movements, a 
mechanism of temporal variations in patch charges is 


m 2 


(•nl 2 )e 2 Dit 

I 6 


(31) 


In Eq. (ED, we have taken into account the dipoles in 
the area ~ nl 2 where n is the dipole concentration equal 
to that of the surface ions. This leads to the estimate, 


6E I Dit 
~e ~ y nl 2 d 2 ' 


(32) 


The relative whisker movement amplitudes A /l are pro¬ 
portional to the above SE/E. As is numerically demon¬ 
strated in Sec. IVIII that ratio is too small to explain 
the observable whisker movements. That smallness is 
due to extremely slow ion diffusion. Even taking ten or¬ 
ders of magnitude higher diffusion coefficient (along grain 
boundaries^) than in Table El retains that smallness. 


C. Non-equilibrium charges 

Some external perturbations, such as light or ionizing 
radiation, create nonequilibrium electron distributions. 
The nonequilibrium electrons drift in the local electric 
fields screening them as illustrated in Fig. [9] These tem¬ 
poral variations in local electric fields can affect whisker 
configurations, which constitutes yet another candidate 
mechanism of whisker movements. It is important that 
such electric currents are much faster than that due to 
the above considered ion diffusion because the electron 
diffusion coefficients are by up to ~ 20 orders of magni¬ 
tude higher. 

Note that the contact with a substrate will addition¬ 
ally enhance the above mechanism because the latter will 
partially absorb either electrons or holes (depending on 
the chemical potentials of a metal film and its substrate) 
thus increasing the electric nonuniformities in the film. 
Therefore, that mechanism is predicted to be substrate 
dependent. 

The problem of photoinduced response of imperfect 
metal surfaces apparently has not been systematically 
addressed. Given a large number of imperfections in the 
metal grain boundaries, oxides, etc., the time scale of 
establishing new steady state distribution of electrons 
can vary by many orders of magnitude between different 
metal recipes and even between different specimens of the 
same recipe. Lacking metal related data one can look 
into the photoinduced recharging data for semiconduc¬ 
tors, especially semiconductor polycrystalline thin films, 
which are ample in the literature for the past decades. In 
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FIG. 9: Energy band diagram for electrons in a nonuniform metal. Ef is the Fermi level separating in equilibrium the filled 
(gray area) and empty states. Left: Equilibrium state. The bottom edge of the partially Elled band is curved due to the random 
fields in the system. The regions of negative and positive curvature correspond respectively to positively and negatively charged 
local regions. Right: Nonequilibrium state under light (or radiation) exposure. The vertical and tilted arrows show respectively 
the processes of photo-excitation and electron drift. The latter processes generate the nonequilibrium population of the electrons 
above the Fermi level thus leveling out the local electric held distribution. 


particular, the phenomena of surface photoconductivity 
and persistent photoconductivity have been studies for 
multiple semiconductor films (see multiple references in 
Refs. I48l - (5g . In may cases, the photoinduced changes in 
charge density were found to be much greater than the 
corresponding equilibrium values; both positive and neg¬ 
ative signals have been observed in different materials. 

The characteristic equilibration times for photoin¬ 
duced charge transfer in semiconductors vary between 
different systems from microseconds to days. Assuming 
similar amplitudes and equilibration times for charges 
photoinduced in imperfect metal surfaces, Eq. pre¬ 
dicts a possibility of ‘spontaneous’ metal whisker move¬ 
ments strongly varying between different specimens and 
light intensity, and ambient conditions. 

A serious problem with the latter explanation is that 
it is difficult to make it more quantitative for the same 
reasons as it is for the case of semiconductors where all 
kinds of imperfections are involved in recharging phe¬ 
nomena. However, that explanation has the advantage 
of offering direct experimental verifications by changing 
the light intensity and temperature while trying to ob¬ 
serve whisker movements. Another attractive feature is 
that it helps to reconcile various seemingly conflicting tes¬ 
timonies, some of which indicate ‘spontaneous’ whisker 
movements, while others do not. 

We shall end this section with one rather arbitrary 
estimate that shows how given plausible parameters of 
the metal surface, the above discussed mechanism can 
predict significant whisker movements. We proceed from 
Eq. (151) with the following modifications: (i) replace the 
ionic diffusion coefficient with the electronic one, Di —>• 
D e , and (ii) multiply the result on the right-hand-side of 
Eq. (1551) with the ratio Sn/n of the photoinduced electron 
concentration change over its average equilibrium value. 
This yields 

SE Sn /ksTnet /00 ^ 

E ~ n Y neEdl ' [ ’ 


Here we have taken into account the Einstein relationship 
between the mobility and diffusivity, p e = ( e/ksT)D e . 
As seen from the numerical estimates in Table IIIIl this 
mechanism appears plausible (note that Sn/n can be 
larger than one in absolute value based on the data for 
photoconductivity in semiconductors). 


VII. NUMERICAL ESTIMATES 


TABLE II: Whisker related geometrical and material param¬ 
eters. 


Parameter 

Value 

Ref. 

density p g/cm 3 

7.4 (7.1)“ 

51 

Young’s modulus Y GPa 

50 (108) 

51 

sound velocity, u s ,10 J cm/s 

2.7 (3.9) 

51 

diameter, 6 d pm 

0.1-20 

2,3,52-54 

length, Z 

10pm-25 mm 

2,3,52-54 

dynamic air viscosity, r] g/cm-s 

2•10“ 4 

55 

kinematic air viscosity, r//pa, cm 2 /s 

0.15 

55 

patch size, do pm 

0.1-10 

7.8 

electric charge density, n e/cm 2 

10 10 - 10 12 

7,8 

near surface field, Eo V/cm, 

<£> 

O 

i—1 

1 

o 

7,8 

ion diffusion coefficient Di, cm 2 /s, 

io- 16 

56 

electron mobility p e , cm 2 /V-s, 

1000 

57 


“The numbers in the top three rows are for Sn; numbers in paren¬ 
theses are for Zn. 

^Whisker diameters and lengths are mutually uncorrelated and 
are described by log-normal probabilistic distributions! 3 ! 52 ! 53 In 
this paper, we use four pairs of values for d and l (Table mill to 
illustrate the scale of variations in their depending quantities. 

Our numerical estimates are based on a set of param¬ 
eters summarized in Table El Note that the difference 
between the most whisker prone metals, Sn and Zn is 
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TABLE III: Characteristic whisker vibration parameters for 
four different geometries. 


Quantity 

d = 1/im 

Z = 1mm 

d = 1/im 

Z = 10mm 

d = 10//m 

Z = 1mm 

d = 10pm 

Z = 10mm 

DaCm/s, 

Eq. m 

5- 10 2 

0.05 

5 • 10® 

5 • 10 3 

Re a 

500 

0.5 

5 • 10® 

5 • 10 4 

0 J b , s" 1 

Eq. © 

3000 

30 

30,000 

300 

A (f>, rad 

Eq. m 

10" 4 

3 • 10" 4 

10"® 

3•10"® 

t Ac/>? S 

Eq. (25]) 

2•10" 4 

2 

2•10" 8 

2•10" 4 

A/Z,Eq. (EH) 
for w = g 

2•10" 3 

1 

2•10" 5 

0.02 

A/Z,Eq. J29l> 

to- 7 

10" 5 

10" 9 

10" 7 

A/Z,Eq. C©), 6 
cos 9 ~ 1 

10" 3 - 10 

0.1 - 10 3 

CO 

1 

O 

1 

u 

1 

O 

i—i 

10 - 10 5 

5E/E , 

Eq. ([32]) c 

10" 9 

io- 10 

10" 9 

10" 10 

5E/E, 

Eq. (©3l) c 
with Sn/n = 1 

0.4 

0.04 

0.4 

0.04 


a Re = Ivapa/v denotes the Reynolds number. 
b The computed ratios A/Z > 1 are formal and must be replaced 
with 1. 

Tor the observation time t = 1 s. 


not very significant and Sn parameters are used every¬ 
where in Table cm 

Table cm contains the derivative parameters quantita¬ 
tively predicting the ability of whiskers to move due to 
the mechanisms considered in this paper. It presents nu¬ 
merical estimates for several quantities defined in equa¬ 
tions that are specified in the table. Such are the char¬ 
acteristic velocity of air flow capable of whisker move¬ 
ments, the Brownian movement whisker fluctuation angle 
and the characteristic time of its accumulation, relative 
whisker deflections due to external vibrations and garden 
hose effect, and relative field fluctuation due to the time 
dependent ion dipoles, and non-equilibrium electrons. 

The relations between these predictions and the 
whisker movement observations were stated in for each 
of the mechanisms in the body of text accompanying the 
equations specified in Table iim 


1) Minute air currents can easily be responsible for the 
movements of thin (d < 1 //m) whiskers, while move¬ 
ments of thick whiskers (d > 10 //in) require much 
stronger, even turbulent (Re > 10 4 ) air flows. 

2) It is much less likely (if possible at all for maybe es¬ 
pecially constricted whiskers) that the Brownian motion 
can be responsible for the observed whisker vibrations. 

3) External mechanical vibrations from power equipment 
etc. can be responsible for whisker movements in a rather 
harsh environment ranging between tracks moving on ru¬ 
ral roads to labs under construction; again, the case of 
severely constricted whiskers make it more likely. 

4) The garden hose instability can never be a source of 
whisker movements. 

5) Electric forces can be responsible for whisker move¬ 
ments if their temporal dependence is determined by the 
electronic (as opposed to the very slow ionic) diffusion. 
That latter possibility can explain the observations of 
‘spontaneous’ whisker movements, and it offers direct ex¬ 
perimental verifications via various dark-light conditions. 

Certain experiments could help to verify the above 
conclusions. In particular, observing whisker dynamics 
vs. (a) light intensity, (b) ionizing radiation, (c) air ion¬ 
izer/deionizer action, and (d) ambient viscosity ranging 
from that of transparent liquids to vacuum conditions / 58 

As another useful outcome, the above results open a 
venue to more quantitative discussions of whisker prop¬ 
erties. For example, they show how to predict whisker 
behavior on accelerated platforms (say, in auto-industry), 
or how to make simple quantitative estimates of whicker 
properties observing their movements under air flows. 

At the end, the question of well enough documented 
observations of ‘spontaneous’ whisker movements calls 
upon further attention, as opposed to the observations of 
air flow effects that have found direct quantitative con¬ 
firmation in the above. The possibility of modifying the 
whisker growing surfaces by shining light or ionizing ra¬ 
diation may have practical applications. 
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